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I. On the Application of the Kinetic Theory to Dense Gases. 

By S. H. Burbury, F.R.S. 
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The motion of a great number of elastic spheres, when their aggregate volume 
does not bear an evanescent ratio to the containing space, has received little attention 
from writers on the kinetic theory. In what respect, beyond the shortening of the 
mean free path, will it differ from that of the rare medium usually discussed 1 I 
think that the answer to this question is that there exists in all systems, dense or 
rare, a tendency for the spheres to move together in masses or streams, and so to 
diminish the mean pressure per unit of area, and the number of collisions per unit of 
volume and time. And this tendency has an appreciable influence on the form of 
the motion as soon as the ratio of the aggregate volume of the spheres to the con¬ 
taining space becomes appreciable. 

If a part of the system, say n spheres, be at any instant contained in a volume V, 
they have energy, T s , of the motion of their common centre of gravity. And they 
have energy, T,., of relative motion. As the spheres increase in diameter, the ratio 
Tr/T, will be found to diminish on avei’age. But the number of collisions per unit 
of volume and time, given T, or T f + T s , depends on TV, and therefore diminishes by 
the diminution of T,.. 

1. Let M be the mass, c the diameter of a sphere, p the number of spheres per 
unit of volume, p the pressure per unit of surface. Also let us now denote by pT,. 
the energy of the motion of the p spheres relative to their common centre of gravity, 
so that T,- is now the mean value per sphere of this energy. 

We have then, as is well known, 

P- lpT r + . . .(|), 

in. which R is the repulsive force, r the distance between a pair of spheres, and the 
summation includes all pairs in unit of volume. We must first evaluate XXKr on 
the assumption that no forces act except during collisions. 

2. Let q be the relative velocity of two spheres. Let 6 be the angle made with q 
by the line of centres at collision, if a collision takes place. The angle 6 may have 
any value from zero to \ir. As the effect of collision the velocity, cos 6, in the 
line of centres is reversed for each sphere. We may assume this reversal to be effected 
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2 MR. S. H. BURBURY ON THE APPLICATION 

by the constant finite force Mg cos O/^clt acting on each sphere in the line of centres 
during the small time 2 dt. 

Let us define this short but finite time 2 dt as the time during which the two 
spheres are in collision, or the duration of a collision. Then, during collision, if dt 
be small enough, the virial of the supposed force is sensibly constant, and we calculate 
its value as follows. Let X, y, v, be the direction cosines of the line of centres referred 
to any axes. The coordinates of the point of contact shall be x, y, z. Then those of 
the centres of the two spheres are, c being the diameter of either sphere, 

x + |Xc, y + \jxc, z + hyc for one sphere, 

and x — |Xc, y — z — fvc for the other. 

The component forces acting at the centre of the first sphere are 
XMg cos 9\2dt, fM.q cos 6/2dt, vM.q cos 6/2dt. 

Those acting at the centre of the second sphere are the same with reversed signs. 

For two spheres colliding with relative velocity q we find that the virial is 
cM.qcos0/2dt at each instant during collision. We have to multiply this by the 
chance that, given two spheres A and B with relative velocity q, they shall be in 
collision at any given instant. 

It is assumed that we are dealing with a space throughout which T,. is constant, 
and therefore the fact that the relative velocity is q, affords no presumption with 
regard to the relative position of the two spheres. About the centre of sphere A 
suppose a spherical surface described with radius c. An element of that surface is 
27 tc 3 cos 0 sin 6 dd. Upon that element of surface form the element of volume 
2ttc 2 cos 6 sin 9 d9 qdt, And form a similar element of volume on the other side of 
the sphere A, that is, using n — 6 for 6. Then the two spheres are at this instant in 
collision if the centre of the sphere B is within either of those elements. 

Let V be the volume in which n spheres are moving with U constant. Then by 
our assumption B is as likely to be in any part of V as in any other. Therefore the 
chance that A and B, having relative velocity q, are in collision, is 

■im 2 cos 6 sin 0 dO qdt 

y 

The average virial for two spheres with relative velocity q is then at each instant 


'w/) 3 {y — v ) 3 -j- (w — w ) 3 j, 


p r ' 2 47 t« z cos 6 sin 6 qdt Mg cos 0 


Y 


2 dt 


M 


3V 


7rc 3 g 3 
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OF THE KINETIC THEORY TO DENSE GASES. 3 

if u, v, iv, u, v', w', be the component velocities of the two spheres. Now since n'\\ 
is the energy of relative motion of n spheres in volume V, 

= MS — — ~ + ^ ~ + ~ ^ , 

r 2 n 

the summation including every pair, and therefoi'e 

ttRr = f 7 tc 3 y 2nT r ..(2), 

= f ttc s p . ‘In’Yy, since — = p. 

Let f tt(? p = k. Then 

22Rr = k . 2nT r for n spheres in Y, 

22Rr = k . 2pT, for p spheres in unit of volume. 

Substituting this value of 22Rr in (1) we obtain 


p = i(l + K)p% .(3). 

It is assumed in these results that we are dealing with a space throughout which 
T, is sensibly constant. 

3. We see then that p is proportional to T } . + /cT r . The analogy between this 
expression and Boltzmann’s T + y, in which y denotes potential energy, suggests 
that the law of distribution of velocities among our spheres should be, instead of 
e“ ; ‘ T as in the rare medium, e~ h(Tr + kT ' > , or rather, since there may be stream motion as 
well as relative motion, e ~ h<T + “ Tr) . 

Let us further develop this analogy. In the Clausian equation 

%pV = nT r + i 2SRr, 

J22Rr = |^ 3 y 2 T,, 

where nT r is the kinetic energy of the motion of n spheres in volume Y relative to 
their common centre of inertia. Hence 



If the n spheres, being initially contained in volume Y 0 , be compressed into 
volume V, and T r be maintained constant during the process, the work done in 
compression is 

w = \'’?dv = {* Blog A + 1 _ A)j T„ 

B 2 
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4 MR. S. H. BURBURY OK THE APPLICATION 

The first term on the right-hand side expresses the work which would be done 
during the process if the spheres were material points, and no collisions took place 
between them. The second term expresses the amount by which this work is 
increased by the spheres having finite diameter c and undergoing collisions. And if 
V 0 be infinite, this additional work is 

m2 

* ttc s -- T„ or § fcnT r . 

V 


We see then that the term hkT,. in the index of e - '“(T + «t,> re p regen ts -work done 
against collisions in compressing the system from an infinite volume to its actual 
volume with constant T,.. It is analogous to the potential y in the usual expression 
c -fc(T + x), It might not be inappropriate to call xT,. the potential of collisions. 

In order to confirm or otherwise the above suggestion, I proceed to consider— 

4. The distribution of energy in a vertical column of gas ‘when in equilibrium in a, 
field of uniform force, the molecules being equal elastic spheres of diameter c. 

Let the column be an infinite cylinder, f, the force, being parallel to its axis. 

Take a. plane perpendicular to the axis as base, and let s be the height of a point 
above that plane. Then we have, with the same notation as before, 


and, as before, 


dp Ids — — M/p 


(4), 


P = 3 (1 + K) pT r . 


Here T,. is the average per sphere of the energy of relati ve motion, which alone is 
concerned in p. But, in our vertical column, assume for the moment that there is no 
stream motion that need be taken into account, and, therefore, we may write T 
instead of T„ and (3) becomes 

p = 3 (! + k) pT . (3a). 

Now k contains p as a factor by (2). But I will now assume that, not T, but 
(1 ~j-. k) T, is independent of s, so that we may write (l + k) T = 3/2 h with li constant. 
On that assumption (3) and (4) give 

p = p 0 e- ftM/s . ..(5), 

where p 0 is the value of p when s = 0. # 

5. Now consider N spheres crossing the plane s = 0 with u for vertical component 
of velocity. Of these, say N — N' reach the plane s = ds without undergoing collision, 
N' will undergo collision before reaching ds. But for every collision by which one of 
the N is struck out, there will be another collision by which another sphere is 

* See Watson’s ‘ Kinetic Theory of Gases,’ 2nd edition, pp. 56-66. 
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OF THE KINETIC THEORY TO DENSE GASES. 5 

substituted for it, the substituted sphere having the same vertical velocity, but not 
quite the same position on average, as the original one had at the instant of its 
collision. 

6. If the collision be direct, i.e., the line of centres coincide with the relative 
velocity, the substituted sphere is advanced in position through a distance equal to c, 
in the direction of the line of centres, and that without any loss of kinetic energy by 
the action of the force f. If the collision be not direct, there is an advance, and we 
proceed to calculate its value. Let A be the centre of a sphere which comes out of 
collision with u for vertical component of velocity, x the centre of the sphere in 
collision with it, H the point of contact. Let l denote the vector line of centres xA, 



and cos (ul) the cosine of the angle between 1 and s. Then ^ c cos (ul) is the 
projection of HA on the vertical, and c cos (ul) is its average value for all the 
collisions in question. 

7. Again, let A' be the centre of a sphere which enters collision with vertical 
velocity u, x' that of the sphere colliding with it, H the point of contact, and let 
V denote the vector A'x'. Then, evidently ^ c cos (ul') = | c cos (ul). 

There will be as many collisions per unit of volume and time of the one class as of 
the other, and the height of the point of contact H, above the base, is on average the 
same for one class as for the other. Therefore, taking the collisions in pairs, one from 
each class, each pair substitutes A for A' as the sphere with vertical velocity u ; and, 
on average, the substituted sphere A is at the instant of its collision above the 
original sphere A' at the instant of its collision by the distance c cos (ul). We have 
next to show that on average of all collisions of N spheres taking place between 
s — 0 and s = ds, c cos (ul) = k ds, where k = |- nc’p. 

8. Let w be the actual velocity of the sphere A as it issues from collision, so that u 
is the vertical component of co, and cos (uw) = uja. Let iff be the velocity of the 
other sphere as it issues from collision with A, and E the angle between <u and iff. 

Also let q be their relative velocity, so that cos (coq) = “—^ e "* - • 

Whatever be the values of &>, if/, and E, 
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6 MR. S. H. BURBURY ON THE APPLICATION 


and 

therefore 


cos (ql) 


* 7 r 

2 cos 9 sin 9 cos 

0 



cos 9 sin 9 d.9 


2 

3 


C COB ~(ul) = 


9c _ 

cos (uq). 

O 


Again whatever be the values of m, *Jj, E, 

cos \uq) — cos (mco) cos (ojq) 

= — cos (ojq), 

therefore 

- 2c U —I 2c U ft) — ylr cos E 

c cos (id) = —-cos (o>q) = —- 1 -• 

d ft) d ft) q 

We have to multiply this expression by the number of collisions which N spheres, 
each having absolute velocity w, undergo in time dt with other spheres having 
absolute velocity xjj ... \p d\]/; the angle between w and i Jj being E ... E -j- c?E, and 
then integrate according to xf> and E. That number is 


N pvc^q dt f (xjj) dxjt J sin E eZE, 

if pf (xfj) dip be the number per unit of volume of spheres whose velocity is iji... xjj -f d\jj, 
so that f f ( xfj) d\jj = 1. 

Therefore the complete average value of c cos (yd), that is the average vertical 
displacement of the substituted spheres, is 


N/37TC 3 dt f # /('!>) f sin E dE q~ cos (ccq), 


or, since dt — ds/u, 


f Npmd dsCdxf,f(xjj) f' £ sin E riE 
Jo Jo 


q COS ( coq ) 


But 


r*« 

Jo 


sin E c^E 


q cos (aq) 


[* , . ft) —-v/rcos E 

| Sin E dK - Y - = 1, 

Jo ® 


also 


ilt — 1. 


so our result is 


f f M # 

Jo 

N §7rc ;i p ds — i<Nds. 


We know (Watson’s ‘Kinetic Theory of Gases,’ 2nd edition, p. 56) that the 
quantity of vertical momentum transferred across the plane s — 0 per unit of area 
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OF THE KINETIC THEORY TO DENSE OASES. 7 


and time is proportional to p, and, therefore, increases in the ratio l : I + K > when 
the spheres, from being material points, acquire diameter c. We now see that this 
is true for each separate class of the vertical momenta. 

9. Corresponding to N spheres with vertical velocity u at the plane s = 0 at the 
beginning of dt, we have at the end of dt N spheres, the same or substituted, whose 
average height is (1 + k) ds or (1 -|- k) u dt. But their loss of kinetic energy due to 
the ascent is on average M fds for each sphere. It follows that the average loss due 


to an ascent ds is, allowing for substitutions, ' . 

1 “J” K 

10. Let then the number per unit of volume of spheres whose energy of vertical 
velocity is ^ id 
Q = M (1 + k) ^ 
have 


i id 


(id -j- did) be at the base, where s = 0, p Q e m dvd in which 
(A). Then the number per unit of volume which at height s 



for energy of their vertical velocity is (remembering (5)) 

p 0 €~ ms dtde~ h{<i ~ Mfs , 

that is pQe~ m dvd (B). The two classes A and B are equally numerous per unit of 
volume, and since, allowing for collisions and substitutions, the loss or gain of energy 

due to the force f in passing from the base to ds or vice versd is Mfdsjl + k, either 
class can by ascending or descending (the proper number of substitutions taking 
place) replace the other. And the assumed law of distribution of velocities is not 
disturbed by the force f, as spheres pass up and down the column. Now, make 
k — 0 in the above reasoning, and we find it is exactly the reasoning from which, in 
the ordinary case of a rare medium, we conclude that T is constant throughout the 
column. The same reasoning leads in the general case to the conclusion that 
(1 + k) T is constant throughout. 

11. The above results are obtained on the hypothesis that no account need be 
taken of stream motion among our spheres. If, however, there be such stream 
motion, we have to suppose that the N spheres crossing the plane s = 0 were 
members of a large group having a vertical velocity U of their common centre of 
gravity, and that the N spheres have vertical velocity u relative to this common 
centre of gravity. Then in time dt the N spheres will have risen on average the 
distance ds by virtue of their relative velocity u, where u = dsjdt, and a distance 
U dt by virtue of the common velocity U. Their loss of kinetic energy by the 
action of the force f in this ascent is 

M/U dt + M.fds/1 + k, 

or if U dt — ds' 

M f (ds + ds /T + k). 
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And now writing T, for T, and T, for the energy of the stream motion whose vertical 
component is U, we find that T + «T r is constant throughout the column, and we 
may now write T -}- kT,. = 3/2 h. 

12. It follows that we cannot express the law of distribution of velocities among 

V/- + + Vj" 

the spheres in the form Ce -,iT , or Ce “ h ~ 2 , as in the rare medium, with 

T = 3/2/i. The law must be Ce -AQ , in which Q is some quadratic function of the 
velocities. Suppose that for n spheres it is 


Q* -— r rp / 1 * -|- + cOcjV.cf T" &c. , 


with corresponding expressions for the components v and w. 

If we find the mean value of Q M by integrating e ” AQ;i for all values of u v « 3 , &c., 
between the limits ± 00 > we find, there being 3 n variables or 3 component velocities 
for each sphere, 

= Sn/2h. 

But 

3%/2/i — n (T + tcT r ), 

therefore, 

(X. = w.(T + ,cT,). 

This result might be considered to justify the assumption that the law of distri¬ 
bution of velocities is in all cases, including the field of no forces formed by making 
f — 0, accurately expressed as follows. The chance that the velocities of n spheres, 
forming a group together, shall be u L . . . u y + du u k c., is proportional to e'-' m<T + ' lTr >, in 
which 

nT = 4 MS {u* -f v 2 + vr) 

nT r — MS {(u — u') % + (v — v'Y -f- (tv — w'Yj/n. 


But we must remember that the whole treatment is based on the consideration of 
a great number of spheres, so that we cannot safely assume the law to hold when n 
is small. Let us, then, consider the subject from yet another point of view. 

13. I have shown elsewhere (‘ Science Progress,’ November, 1894), that in a dense 
medium the velocities of contiguous spheres cannot be independent of one another, 
because there is a presumption that recent collisions of two spheres near to one 
another have been with the same third sphere, and they have, so to speak, inherited 
some common velocity from it. In other words contiguous spheres have been exposed 
to the same environment, and, in the dense medium, environment does not change 
rapidly ; therefore their velocities are not independent. 

If, therefore, the spneres contained in a volume V be n in number, and their 
positions known, the chance that their velocities shall be 
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must be of the form 


u l + du\, &c., tq . . . v x + dv lt &c., 

h («iRi 2 + h% u i l h + a 2 % 2 + <&«•) 


9 


in which the index is a quadratic function of u v u 9 , &c., v v v 2 , &c., and iv 1} iv 9j , &c., 
but contains no products of the form uv, uw, or vw, and the coefficients b l9 , &c., are 
functions of the positions of the spheres. 

14. With regard to the forms of these coefficients, we observe that the quadratic 
function in the index must always be positive, because the chance cannot become 
infinite for infinite values of the variables, u. The condition for this is that the 
determinant 


D = 


2oq, b 12 , btf. 

^12> 2a 3 , & 33 , 




and all its coaxial minors, must be positive, and, therefore, every a positive. 

15. Again b n , not being zero, expresses the fact that u l and u % are not independent. 
But it is also a fact that they are more likely to be of the same than of opposite 
signs. Therefore, b n , and similarly every b, must be negative or zero. 

16. Evidently also the coefficients b must generally diminish in absolute magni¬ 
tude as the distance between the spheres to which they relate increases, and 
must become inappreciable at some distance, small compared with the dimensions of 
our system, but possibly large compared with the diameter of a sphere. The b’s 
must be functions of the positions of the n spheres within V having this property. 

Again we may consider either a volume Y containing n spheres, or a smaller 

7b 1 

volume-- V containing n — 1 out of the n spheres, that is, all the n spheres 

7b 

except one, which one belongs to the outer layer. If the velocities of that one 
be denoted by u„, v„, w n , and if Q„ be the quadratic function for the n spheres, Q„_ 
for the n — 1 spheres, we must have 

+ oo 

£ -«q„_x _ |j|€ _/ ‘+ du, i , dv u , dw n . 

— ao 1 

If we actually perform the integrations indicated for one variable, we get with 

Q, t = apq 3 + + &c. 

Q«_i a ]tq 3 + b + co + &c., 

in wdiich 

2a' x = 2 a l — b v “\2a„, 2 a' % — 2% — b 9 ,?j2a m 

b\ 9 — b l9 — & 1) Ao„/2a„ &c. 

MDCCCXCVI.—A. c 
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10 MR. S. H. BURBURY ON THE APPLICATION 

As the number of spheres included in our group diminishes, the coefficients a 
diminish, and, since every b is negative, the b’s increase in absolute magnitude. 

On the other hand, as n increases the a’s increase, and every 5 s , if changed at all, 
diminishes. Therefore, as u increases the function Q„ tends to a limiting form. 
But that limiting form must be »T if k = 0, because in that case we know the 
law to be e -7lT . We must then have ultimately Q„ — n (T + k6), where 0 is a 
quadratic function of the velocities. And we may now assume that for a sufficiently 
great number of spheres comprised in a group, throughout which k is sensibly 

1 

constant, the law is expressed by the function Ce "»< T + * T r> j in which C = 

D being the determinant of the coefficients of the quadratic function T + kT,.. 

We have thus obtained certain conditions which the coefficients a x , b 12 , Sec., must 
satisfy. Another condition is that the assumed law of distribution of velocities 
expressed by the function 

(c^uf + -b &c.) 

shall not be disturbed by a collision taking place between any two of the n spheres, 
which collision changes the velocities, but not the positions, of the two spheres in 
question. To find the values of a u b n , See., to satisfy this condition presents consider¬ 
able difficulty. It is not, however, at present necessary to solve the problem in that 
form, as will be seen later. For we have only to consider the positions of the n spheres 
as unknown, and we obtain a solution sufficient for our purpose. If, namely, it be 
given that there are at any instant n spheres within a spherical space S, but nothing 
is known of their positions within S, we have only to assume that the chance of their 
having at that instant velocities u 1 . . . u x + du x , &c., is 

e~ 7iQ du x . . . dw n , 

with 

Q — at ( u 2 + v 3 + tk 2 ) + bhh (tiu' + vv + ww') 

containing only one coefficient a and one 6, and we shall find that all necessary 
conditions are satisfied, including the condition that the assumed distribution shall be 
unaffected by collisions. 

For let x } , y v z v aj a , y. 2> z 2 be the component velocities of two spheres before collision. 
A collision between the two converts these components into x' x , y\, z\, x' 2 , y' 2 , z\ in 
the following manner. Let X, /x. v be the direction cosines of the line of centres at 
collision. 

The velocities of the two spheres resolved in the line of centres are, before 
collision, 

\x x + /xjq + vz x and kx. 2 + yy 2 + rz 2 


respectively. And we have 
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X \ — x \ — A (^1 + Wy 1 + vZ i) + ^ (^ X -2 + Ms + vz 2 ), 

y\ — V\ ~~ p (^ x i + yy 14- vz x ) 4- y (Aa% 4- yy% 4- v h), 

&c. 

That is, 

x\ = (1 — A 2 ) x 1 — kyy x — \vz x + A 2 » a + A/xy 2 -j- Ar% 

y'i — — + 0 — / xS ) ?/i — + V *2 + 

z\ — — \vx x — \mj x + (1 v 1 ) z x -\- \vx 2 -f jxvy % + v\ 

x' 2 — A 2 aq -f- \yy x A vz x + (1 — A 2 ) x 2 — A yy. 2 — A vz 2 
y\ — V *1 + M 2 2/i + yvz x — A/xjc 2 + (1 — /X 2 ) y 2 — y,vz z 
z 2 — \vx x 4~ l xv lh + — kvx 2 — /xry 3 4” (1 — r 2 ) Ba¬ 

call these equations A. 

But inasmuch as the motion might take place the reverse way with the same 
values of A, /x, v, it must by the same reasoning be true that 

x x = (1 — A' 2 ) x' L ■— \[xy\ — A vz\ 4- A 2 ,^4 4* A/xx/' 2 4- Aaz' 2 , &c.A', 

which are the same as equations A with the accents interchanged between the right 
and left-hand members. If we solve either set of equations we get the other set. 

17. Now it is given that at this instant the chance of the spheres having velocities 
4- dx x , &c., is Ce~ 7iQ dx x dx 2 ... in which 

Q = ax i 4" bx x x % + ax 2 4- ay x 4" by x y 2 + a Vi + az \ + bz x z 2 4~ ct2 3 2 

4- bx x (x 3 4" x 4 + &C.) 4~ &c. 

4- cix./ J 4~ &c. 4~ terms containing squares and products of the 
velocities x s . . . z n . 

The chance that after this collision the velocities of the n spheres shall have the 
values x\, y\, z\, x' 2i y\, z' 2 , and x S) y 3 , &c., is found by substituting for x x , y x . z x , 
x 2 , y 2 , z. 2> in the index AQ, their values in terms of x\, y\, z\, x' 2> y\, z 2 , as given by 
equations A'. Effecting this substitution we find that the coefficient of x'^ is 

a . {(1 - A 2 ) 2 4- AV 2 4- AV 4- A 4 4- A 9 > 2 4- AV 2 } 

4 -b. {A 2 (l — A 2 ) — Ay — AV] 

that is a. Similarly, we find that the coefficient of x\x 2 in the new index is b, and 
the coefficient of every product x\x r in the new index is the same as that of x x x r in 
the original index. The new index is then the same function of x\ , y\, z\, x\, y' 2 , z\, 
that the original one was of x x , y x , z x , x 2 , y 2 , z 2 . The assumed, distribution of velocities 
is therefore not affected by any one, and therefore not by any number of collisions. 
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12 MR. S. H. BURBURY ON THE APPLICATION 

This, it will be remembered, is on the assumption that all the coefficients a are the 
same, and all the coefficients b are the same, so far as they have to do with the 
colliding spheres. 

18. Now 


M 


% (T + kT t ) 

n — 1 \ uf 


1 + 


K ) 


2 


k ic , n — 1 \ u 2 , s 

n 1 2 n 1 J 1 \ n ) 2 1 

n — 1 


at (u~ + v 2 + w 3 ) + btt. (uu! + vv + low'), if 2a = 1 +- K and h 


IC 

n 


and if we use this for Q« we find that all necessary conditions are satisfied, including’ 
the condition for permanence notwithstanding collisions. 

19. If, however, the coefficients a, b were not the same for both colliding spheres, but 
the form were a x x x + b^x^ fi- %a: 2 2 -fi b^x^c,^- b 2 ^x 2 x 3 -\- & c -> then we should find that 
the coefficient of x in the new index is (1 — X 2 ) a 1 + X 2 » 3 ; the coefficient of *' 2 2 is 
(1 - X 2 ) a z + Xht}; the coefficient of x\x\ is as before, b n ; but that of x\x 3 is 
(1 — X 2 ) b is + X 3 &23 > and that of x\x 3 is (1 — X 2 ) b 2S + X 2 & 13 . 

The assumed law of distribution cannot in this case be unaffected by collisions, 
unless (l) all the a coefficients are the same; (2) if the velocities of the two colliding 
spheres be u u v v u\ and u^v^Wa, and those of any third sphere be iu, v 3 , iv 3 , then & ;3 =& 23 , 
that is, the b coefficients must be such functions of the position, that if spheres 1 and. 
2 are close together 6 13 = b 2 s , &c. 

To return to the case of the positions of the spheres being unknown. 

20. If we form the determinant of the system 





K 


Vv 




K 

n 

K 

n 


we find 


Hence w’e find 


D = (i + k)'* — (l + k )” -1 
Dn = (1 ~r «) n ~ l ~~ n ~ L " (! + «)‘ 


:r 3 


•-ft p 

tlJ = &C. 


I) n _ 1 1 . n + k 

D = AM ’ n ‘ TJlc ’ 


therefore 
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n + k 
1 -f* tC 


u i u -i = > &c. 


13 


[These results are easily obtained by considering the general determinant of n 3 
constituents 

2 a b b b . | 

D = b 2a b 


in which all the axial constituents are 2a, and all the non-axial constituents are b . 
It will be found for n = 2, w = 3, and thence by induction, that 

D = (2a ~~ ft)* + nb (2a — 6) n “h 


Replacing 2a by 1 + n —■ 1/e/n, and 6 by — k[u, we get the result above stated.] 

Also and w 2, have the same value as v?. And the whole kinetic energy of the 

, m . ,«• 3 nv? m 3 n -F tc 

n spheres, or ni, is, on average, 1VI, or n 1 = — Y+le * 

Again, 


k rn w —• 1 

7"™ n ±r = K “X—. 

AM 2n 


— UAl 9 — — 
n 1 * n 


u Y u 3 — &c , 


+ K 


n — 1 


2n 


Ucy 


71 


UM S — &C., 


— &C., 

with similar expressions for the v’s and w 9 &. But 
n — 1 


1-1 + K 




7 — — — &c. = (i + k ——b bi 

n ' \ % J D 


J c Bi 3 

~n D 


&c. = Sn/2KNL 


by the properties of determinants. 
Therefore 


and 


KfhVr 


3 n 
2Ti 


on D n 
21b D 


3 n 
2 h 


3 n + tc 
2 h 1 -f tc 


nT r 


3 n — 1 


2h 1 *4~ /c 


3 w 

2 lb 1 tc 


nT s = nT — wT r 


3 In tc 
2 h \1 -|- tc 


n — 1 
1 ~f- tc 



and therefore 
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and therefore T.,/T — (1 + k) j{n + k), which increases as k increases, that is, 
cceteris 'paribus, as the diameter c increases. 

Again, the pressure per unit of area on a plane moving' with the stream, and 
therefore the mean pressure, is p — f (1 -f- k) pT„ which is independent of the 
diameter c. 

Now, the spheres being material points with T r = 3/2 h, the mean pressure is 
| pT, ; that is, as we may write it, f-(l + k) pT r , because in this case k — 0. As 
the spheres increase in diameter with (1 + k) T,. unaltered, the mean pressure per 
unit of area remains unaltered. In other words, it is exactly as much diminished by 
the conversion of part of the energy of relative motion of contiguous spheres into 
energy of stream motion as it is increased by the introduction of the term SSKr as 
the spheres acquire diameter c. 

Comparing the actual value of p with what it would be if, with the same total 
kinetic energy, the spheres were material points, we see that it is diminished in the 
proportion 1 + ^:1. 

The number of collisions per unit of volume and time is proportional to c 2 T,.; that 
is, to k*/(1 + «). It is less in the proportion 1/(1 + k) than it would be if, with 
the same total kinetic energy, the spheres had velocities independent of one 
another. 

21. From the fact that p is independent of k, it follows that local variations of 
density, i.e., of k, do not involve the expenditure of any work on the whole, and 
therefore such variations may and will come into being. 

22. Boltzmann’s minimum function continues to diminish by collisions, finally 
attaining its minimum constant value when the distribution of velocities defined by 
our assumption in 13 is established, but its actual value when minimum differs by 
^ log D from what it would be if the spheres were material points. 


